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We present that a linear response theory in the continuum can be easily formulated 
with Absorbing Boundary Condition (ABC). The theory is capable of describing continuum 
spectra and dynamical correlations. Application of the ABC does not require the spherical 
symmetry and the method is suitable for mesh representation in the real coordinate space. 
Isovector giant dipole resonances in beryllium isotopes are studied with the time-dependent 
Hartree-Fock with the Skyrme force in a three-dimensional mesh space with the ABC. 
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mation (CRPA). 1 ) The trick to treat the continuum was use of the single-particle 
Green's function with the Outgoing Boundary Condition (OBC). This relies on ex- 
plicit calculations of Green's functions with the OBC in the one-dimensional (radial) 
coordinate. We have developed a method to calculate the continuum response in the 
three-dimensional (3D) real space by solving equations iteratively with the Green's 
functions. 2 )' 3 ^ The OBC is explicitly taken into account with this method. However, 
it is costly in computer resources. We have also been investigating a possibility to 
approximate the OBC by Absorbing Boundary Condition (ABC). 2 )" 4 -* An advan- 
tage of the ABC is its simplicity and flexibility. For instance, it is extremely difficult 
to treat the OBC in the real-time calculation of the time-dependent wave functions, 
while it is straightforward to implement the real-time method with the ABC. 2 ) 

The Skyrme Hartree-Fock (HF) theory 5 ) has been extensively applied to study 
of ground-state properties of relatively heavy nuclei. To investigate their excited 
states, a straightforward extension is the linear response calculation based on the 
HF ground state. The CRPA combined with the Skyrme Hartree-Fock (HF) theory 
has been extensively utilized to study giant resonances in spherical nuclei. 6 ) " 8 ) So 
far, only spherical nuclei were studied because of the limitation of the conventional 
CRPA. We would like to show that the linear response calculation with ABC is a 
feasible extension of the spherical CRPA to the deformed one. 
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§1. Introduction 



J> \ The drip-line nuclei are finite fermion systems whose separation energy is nearly 

zero. For such weakly bound systems, the continuum should be properly taken into 
account in description of their structures and reactions. In the linear response regime, 
the inclusion of the single-particle continuum for particle-hole (p-h) excitations has 
been achieved by using a method, so called, Continuum Random-Phase Approxi- 
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§2. Skyrme Time-Dependent Hartree-Fock (TDHF) with ABC 

Since the Hamiltonian in the Skyrme HF theory is almost diagonal in coordi- 
nate representation, a mesh representation in the 3D coordinate space provides an 
economical description. The real-time calculation of the Skyrme TDHF were carried 
out for studies of heavy-ion reactions. 9 ) The main issue is then how to incorporate 
the proper boundary condition in the 3D uniform grid representation. In this paper, 
we adopt a method of ABC combined with the real-time TDHF. 2 ) 

As is well known, the Green's function with OBC (for E > 0) can be simply 
written as 

G^Hr,r>;E) = (r\ E _ 1 H + i7] \r>), (2-1) 

where rj is a positive infinitesimal. In numerical calculations, we take a spherical 
box with radius R as a model space. Since the infinitesimal r\ cannot be treated 
directly in actual numerical calculation, we need to explicitly construct the outgoing 
solutions at the boundary (r > R). This is a difficult task when the system does 
not possess a spherical symmetry. The idea of ABC is to replace the rj by a finite 
coordinate-dependent quantity, (r|^|r') = fj(r)S(r — r'), and imposes the vanishing 
boundary condition: 

G ^c) {r y. E) = (r | 1 | r '). (2-2) 
E — H + iT] 

Here, fj(r) is taken to be positive far outside the system (R < r < R + Ar) 
and zero elsewhere. The Green's function has a vanishing boundary condition, 
G^ ABC \r,r'; E) = for r = R + Ar or r' = R + Ar. This is equivalent to bound-state 
Green's function for the Hamiltonian with a complex absorbing potential, H — ifj(r). 
Now let us discuss properties of G^ ABC \E). 

For E < 0, the Green's function with ABC has a pole at eigenenergies E = E n , 
because the bound eigenstate is localized, ipn{ r ) ~ at r > R. For E > 0, the 
Green's function has an outgoing asymptotic behavior. This can be seen as follows: 
In the outer region of our model space (r > R), the complex potential —if] produces 
complex wave numbers, k + i~f where 7 > 0. Thus, the outgoing waves exponentially 
damp while incoming waves diverge. The vanishing boundary condition selects only 
the outgoing solutions. In the inner model space (r < R), as far as the reflection 
caused by the potential, —ifj(r), is neglig ible, G( AB °)(£) is identical to G^ 0BC \E). 
Therefore, the complex boundary potential must be strong enough to absorb the 
whole outgoing wave and simultaneously gentle enough to minimize the reflection. 
We adopt an absorbing potential of linear dependence on the radial coordinate, 

~( \ - / °' iorr< R, 

V[r) ~ {VO^, for R < r < R + Ar. ^' 6 > 

In order to minimize the reflection, the height r/o and width Ar should satisfy a 
condition 10 ^' 2 ) 
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<7] < -^-ArV8r~^E 3/2 . (2-4) 
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An advantage of the ABC is its applicability to the real-time TDHF calculations. 
We use the Skyrme TDHF method to investigate giant resonances in the continuum. 
The initial state at t = is prepared by boosting the HF ground state, |#b)> with 
an instantaneous external field, V e ^t{r,t) = v e ^t{T)5{t). Then, the time evolution of 
the TDHF sate, \\P(t)), in a 3D real space is calculated according to a prescription 
given in Ref. 11). Time step is taken as At = 0.001 MeV -1 and iV-times iteration 
of time evolution is done for ^(t)} up to the total period T = NAt. After the 
time evolution, we carry out the Fourier transform of expectation values of a one- 
body operator D, (# r (i)|l)|<P , (t)). The imaginary part of this quantity gives the 
corresponding excitation strength, S(D;E) = J2 n \(n\D\0)\ 2 5(E — E n ). See Ref. 2) 
for details of the numerical calculation. 



§3. Numerical results and discussions 



3.1. Continuum effects: an illustration 

In many cases, the RPA response 
functions have been calculated utilizing 
the L 2 basis set, such as harmonic os- 
cillator basis. A similar L 2 -type calcu- 
lation is possible using the coordinate 
space of finite box with the vanishing 
boundary condition. We call this "Box 
Boundary Condition (BBC)". The low- 
energy RPA response below the separa- 
tion energy are well described with the 
BBC. Then, in order to obtain a rea- 
sonable BBC description of the response 
in a energy region above the separation 
energy, what size of model space is re- 
quired? In other words, how large the 
model space should be, to simulate the 
continuum wave functions? 

Assuming the BBC of box radius R, 
an escaping particle of velocity v (en- 
ergy E) reflects at the boundary wall at 
r = R. This reflection gives a spurious 
contribution to the continuum response 
of the system and eventually produces a 
discrete spectrum. Thus, the maximum 
period of physical time evolution is 

T = R/v = R^MjlE. (3-1) 

The uncertainty principle tells that the 
total period T is related to an energy 
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Fig. 1. Strength functions of giant monopole 
resonances in 16 O calculated with the BKN 
interaction. A smoothing parameter, F = 
0.5 MeV, is added to the energy as E + 
ir. (a) Results of CRPA using a radial 
coordinate up to R = 10 fm. (b) RPA with 
the BBC of R = 10 fm. (c) The same as 
(b) except for R — 100 fm. (d) The same 
as (b) except for R = 200 fm. 
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resolution AE after the Fourier transform: 



AE = — 
T 




(3-2) 



Here, M and E are the mass and the kinetic energy of escaping particles, respectively. 
Therefore, for instance, if we require AE = 1 MeV for the strength function at 
E = 10 MeV above the threshold, one has to handle the box of R w 200 fm. Dealing 
with this large space is difficult even for a spherical (ID) case and is almost impossible 
for a deformed (3D) case. 

Using a simplified Skyrme force, the BKN interaction, 12 ) monopole response in 
16 is calculated both with and without the continuum. For this case, owing to the 
simplicity of the BKN interaction and the spherical symmetry of the HF state, we 
may carry out the RPA calculation with the BBC in a large coordinate space. We 
also perform the CRPA calculation. The OBC calculation is done for a spherical 
box of R = 10 fm, the result of which is compared with the BBC results with three 
different boxes of R = 10, 100, and 200 fm (Fig. 1). The figure shows that we need 
utilize a coordinate space of R > 100 fm to obtain a sensible result. On the other 
hand, if we use the complex potential of BBC, Eq. (2-3), a space of Ar = 10 fm in 
addition to R = 10 fm is enough to obtain an identical result to the OBC. 4 ) 

3.2. Giant dipole resonances in beryllium isotopes 

Isovector dipole states in even beryllium (Be) isotopes ( 8_14 Be) are investigated 
using the real-time-ABC method in Sect. 2. We adopt the SIII parameters for the 
Skyrme interaction including all time-even densities and most of time-odd densities. 
!3) Properties of the HF ground states of these isotopes were investigated in Ref . 14) . 
The ground-state binding energies are well reproduced with the SIII and prolate 
shapes of the density distributions are obtained for 8,10,14 Be. We now investigate the 
El linear response on these ground states. An initial state for the real-time TDHF 
calculation is prepared by boosting the ground state with an El operator with a 
recoil charge, t> cx t = (Z/A)r n — {N/A)r p . The total period of time evolution is taken 
as T = 4 MeV -1 . This corresponds to about an energy resolution of 1.6 MeV. The 
model space we take is a 3D lattice coordinate space of R = 8 fm with a square mesh 
of Ax = Ay = Az = 1 fm. For the ABC potential, Eq. (2-3), we take r? = 14 MeV 
and Ar = 10 fm. 

In Fig. 2, the obtained El strength functions are presented. Since 8 > 10 > 14 Be h ave 
prolate shapes in the ground states, the El strengths parallel to the symmetry (z) 
axis and perpendicular (x and y) axes are displayed separately. Large deformation 
splitting over 10 MeV is obtained in these nuclei. This is because the ground states 
have large deformations which almost correspond to super deformation in heavy nu- 
clei. The ground state of 12 Be is nearly spherical and we show the El response 
summed with respect to three directions in Fig. 2 (c). 

Appearance of low-energy (pigmy) strengths (E < 10 MeV) is seen in 12 Be 
and even enhanced in 14 Be. These strengths disappear if we get rid of the ABC 
potential, which suggests that the continuum is essential to produce the low-energy 
El strengths near threshold. It is also worth noting that the real part of the response 
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Fig. 2. El strength functions of giant dipole resonances in even- A Be isotopes, calculated by the 
TDHF with the SIII force using the ABC potential: R = 8 fm, Ar = 10 fm, r) = 14 MeV. 
El strengths parallel (perpendicular) to the symmetry axis are shown by dashed (dot-dashed) 
lines. 



function changes its sign at the position of the low-energy peak. This may suggest 
a real resonance character. 

§4. Conclusions 

In this paper, we present the method of ABC which is able to describe the 
3D continuum effects in a feasible manner and in good accuracy. We take El giant 
resonances in beryllium isotopes and apply the real-time TDHF- ABC method. Large 
deformation splittings in 8 ' 10 ' 14 Be and low-energy resonances in 12,14 Be are obtained. 

Although we have shown only the application to the response calculations in 
this paper, the ABC provides a general and powerful method to treat the OBC. 
Another application, a calculation of breakup processes of a halo nucleus, is also 
under progress. 15 -* We are encouraged by the results we have obtained and will 
study structures and reactions of unstable nuclei using the ABC. 
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